For Reference 


NOT TO BE TAKEN FROM THIS ROOM 





For Reference 





NOT TO BE TAKEN FROM THIS ROOM 


si 1IBRIS 
WRRSITAIS 
atk THESIS 






\ I 
S| = = 
aa | sam 
| 








DBS 3803 

The University of Albert 
Printing Department 
Edmonton, Alberta 














f 


Tan WetvaNatty or ALanetA 


STUDIES ON ACKLERRAR Other t-1ATA 
SYSTEMS Wo NONHEREO IS fPLiL “ceerricths 





; | a 


(¢) Mi) QMARED RIA) ®.n0- 


a ae 
SUEMeTTED cm Mee « 


IN PARTIAL SEZice wh 





Digitized by the Internet Archive 
in 2022 with funding trom 
University of Alberta Libraries 


https://archive.org/details/Rashed1970 


THE UNIVERSITY OF ALBERTA 


STUDIES ON NONLINEAR DISCRETE-DATA 


SYSTEMS WITH NON-ZERO INITIAL CONDITIONS 


by 


AHMED MOHAMED HELMY RASHED 





A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 


OF MASTER OF SCIENCE 


DEPARTMENT OF ELECTRICAL ENGINEERING 


EDMONTON, ALBERTA 


SPRING, 1970 


“y 


ry 


TORNTA “UW 





| 
; : 
mt 9 _ 
: - 
. 
| wed 

, : 
7 

~, 

7 a 

| : ae OP, "Ly Dal 





UNIVERSITY OF ALBERTA 


FACULTY OF GRADUATE STUDIES 


The undersigned certify that they have read, and 
recommend to the Faculty of Graduate Studies for acceptance, 
a thesis entitled "Studies on Nonlinear Discrete-Data Systems 
with Non-Zero Initial Conditions" submitted by Ahmed M.H. Rashed 
in partial fulfilment of the requirements for the degree of 


Master of Science. 


a & YP 





. { w=; io &!? Beeeraheg 


ions } [ . i Waee hy ij i Te 3 6° 





(i) 


ABSTRACT 


Nonlinear discrete-data systems subject to non-zero initial 
conditions are considered. Two types of systems are investigated; the 
first has no hold device, while the second includes a zero order hold. 
In both systems the nonlinear element is located in the feedback path. 
Linear discrete-data systems with non-zero initial conditions are 
considered first and then it is shown how one, similarly, can treat 
nonlinear discrete-data systems with non-zero initial conditions. In 
the latter case system response takes the form of the discrete Volterra 
series in which non-zero initial conditions are included. The series 
is obtained by solving the system equation by iteration. Uniqueness 
as well as convergence of the series, and consequently bounded-input 
bounded-output stability of the nonlinear system, are investigated by 
use of the contraction mapping principle. By letting the system in- 
put be zero, asymptotic stability of the system can be studied. Finally, 
the intersampling response of nonlinear discrete-data systems is obtain- 
ed in the form of an infinite series, and the stability of the system 


between sampling instants is discussed. 
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CHAPTER I 


INTRODUCTION 


A control system in which the signal at one or more points 
is in the form of pulses occurring at discrete instants of time is 
called a sampled-data or a discrete-data system. If the time in- 
terval between every two consecutive pulses is a constant, the signal 
is said to be uniformly sampled. A typical sampled-data system 


utilizing uniform sampling is shown in Fig. 1.1. 
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Discrete-data systems have drawn the attention of many in- 
vestigators since the Second World War. The main reason is the rapid- 


ly growing use of digital equipment, mainly digital computers, in 
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modern scientific and technological applications. Examples are: 
computer-controlled machine tools, time sharing of transmission links, 
ground-controlled guided missiles, etc. 

An important field of investigation into the behaviour of 


these systems is system stability. However, much of the work done 


OFPthe mest ’sienificant*work inthis context is that of Meyeieiae a who 


extended the results of V.M. Popov for nonlinear continuous-data systems 


to discrete-data systems. Jury and Retwres further extended the work 


of Tsypkin, and by constraining the slope of the nonlinearity, they 
were able to obtain less conservative results. 

Another important concept of system stability is the concept 
of bounded-input bounded-output (BIBO) stability. Here a system is 
said to be stable in the BIBO sense if, and only if, every bounded 
input produces a bounded output. For linear systems the equivalence 
between system stability in the BIBO sense and in the Liapunov sense 

: 133 4 ; 
has been established . On the other hand, although it has been 
shown, for a certain class of nonlinear systems, that conditions suf- 
ficient for establishing asymptotic stability also establish BIBO 
ears el ee ead Re, ; F ; F 
Stability , there is, in general, no simple equivalence between 
1 
these two concepts and counter examples have been constructed af 

Volterra series provides a powerful tool for finding system 

19=24 


response and system stability of nonlinear systems we nesconvers= 


ence of the series implies that the system is stable. 


* Numbers placed above the line of text refer to references. 
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in this area so far discusses only the problem of asymptotic janie a 
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It is the objective of this thesis to use Volterra series 
to study ayant discrete-data systems with non-zero initial 
conditions. The response of the system at sampling instants, as well 
as in between sampling instants, is represented by a discrete Volterra- 
type series. BIBO stability as well as asymptotic stability are then 
studied through the convergence of vs resulting series utilizing the 
Banach fixed point Siena 

A brief review of the important results obtained so far in 
the area of system stability is given in chapter II. In chapter III 
linear discrete-data systems are considered. Non-zero initial condi- 
tions are included in a general expression for system response. Both 
BIBO and asymptotic stability are discussed using functional represent- 
ation of the system. The functional notation and concepts involved 
are introduced in that chapter. In chapter IV the development of 
chapter LII is extended to nonlinear discrete-data systems. The res- 
ponse of the system is obtained in the form of a discrete Volterra-type 
series by a process of successive approximations. Chapter V deals with 
the response of the nonlinear discrete-data systems between sampling 
instants. 

A brief review of Z-transform theory, which is the main tool 
“used to associate initial conditions with system response, is given in 


the appendix. 
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CHAPTER II 


REVIEW OF DISCRETE-DATA SYSTEM STABILITY 


Zoe LN Ee GOduCctLon 

In considering the stability of a given control system we 
often consider both asymptotic and bounded-input bounded-output (BIBO) 
stability. An undriven system is said to be asymptotically stable if, 
given some set of initial conditions, its output tends to zero as time 
tends to infinity. On the other hand a driven system is said to be 
BIBO stable if its output remains bounded for each bounded input. 
However, while it is true that for linear systems asymptotic stability 
implies BIBO stability, this is not, in general, true for nonlinear 
ausuaende 

The purpose of this chapter is to review briefly the import- 
ant results obtained in the field of stability of nonlinear discrete- 
data systems. The first section discusses asymptotic stability results 
and the second BIBO stability results. 
2.2 Asymptotic Stability of Nonlinear Discrete-Data Systems 


The first extensive work in the field of asymptotic stability 


>“ and Kalman 


of nonlinear discrete-data systems is the work of Hahn 
* 

and Rowan . They extended the ideas of the second method of 

Liapunov to apply to systems of difference equations. The results 


they obtained governing the choice of Liapunov functions suitable for 


different types of stability of the null solution are analogous to 





* See also Ref. 4. 
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those obtained for continuous-data systems. Although this method is 
very general, its application requires both intuition and experience 
in the choice of suitable functions. 

The same problem was treated by the Russian scientist 
Ya.Z. Tsypkin who has shown that it can be solved more easily and 
effectively by working in the frequency Aone He extended the 
ideas of the Rumanian scientist V.M. Popov for the determination of 
the stability in-the-large of continuous autonomous systems to dis- 
crete-data systems. This method yields results which are intimately 
connected with the usual notion of frequency response characteristics 
and, more importantly, gives general sufficient conditions for stabil- 
ity which are applicable to systems of arbitrary order. 

The system considered by Tsypkin consists of a nonlinear 
element located in the forward path of an otherwise linear feedback 
system, Fig. 2.1(a). The nonlinear characteristics $(x(n)) is assumed 
to be memoryless and confined to a specified gain sector [0, k), 

Fig, 2.1(b). Mathematically the nonlinearity is assumed to satisfy 
the conditions: 


eee = 0 
C21) 


_ x(n) 
tee ee Go 





Figs 2.1) (a) 
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Fig. 2.1 Nonlinear discrete-data system considered by Tsypkin 
(a) system configuration 
(>) nonlinearity characteristics 


The results obtained by Tsypkin are similar in many respects 
to the Nyquist criterion for the stability of linear discrete-data 
systems. However, they give only sufficient conditions for system 
Stability. The system of Fig. 2.1 is asymptotically stable in the 
large if the following inequality is satisfied: 


Ree! oe +i 


= 10 CZ) 
where 


c(t) = G(z)|_ jul (2.3) 


T is the sampling period. 
Geometrically, inequality (2.2) means that the frequency 
characteristics mee of the linear part must be located to the 
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plane at the point 5, 


Os 

Tsypkin further extended the class of systems for which in- 
equality (2.2) gives sufficient conditions for Stability. Time-varying 
nonlinearities and systems with neutral or unstable linear part are 


included. In this case the nonlinearity, instead of satisfying condi- 


Eton. Im) baseto satasty the condition: 


- < ¢@),n) 


ep) <kte | (Qn) 


where © is a real number stuch that 


os PhaaG Cz) 
ce 1+eG(z) Nee 


papctap le: and (es) takes the place of eed sey MOLES) 5 

Tsypkin also noted that for some cases condition (2.2) is 
not only sufficient but also necessary for asymptotic stability. 
These cases arise when the maximum negative real part of dee, or 
ye iovattained at the abscissa. This is becausescondition (2.2) 
in this case coincides with the necessary and sufficient conditions of 
Stability of the linearized system. 

tsypkin, inyaolater ae formulated his results into the 
familiar cirele criterion. 

The work of Tsypkin was further extended by Jury and Tecgaaae 


They considered the same class of systems considered by Tsypkin with 


the addition that the nonlinearity is to satisfy the extra condition: 


COC) eee tee 
| GD | k' < G20) 


or more generally 


do(x(m)) . 
ES Se ERT ERY it 1S (2.7) 
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where k' is a positive real number; and ky and ky are any real numbers. 

The approach used by Jury and Lee is also a frequency domain 
approach and, like Tsypkin's approach, depends ae the frequency 
characteristics of the linear plant which is assumed to have all of 
its poles in the left half of the s-plane except for the possibility 
of one pole at the origin. 

Like Tsypkin's results, the results arrived at by Jury and 
Lee give only sufficient conditions for asymptotic stability in-the- 
large of the nonlinear discrete-data system belonging to the class 
considered. These results are given in the form of two inequalities 
each of which applies to certain types of systems. 

For a nonlinear discrete-data system whose nonlinearity 
Satisfies conditions (2:1) and (2.6) the Satisfaction of the follow- 
ing inequality on the unit circle |z| =IPisee Guiticlene cOnai ton 


for asymptotic stability in-the-large of the system: 
! 
Re (6(z) [1+q(z-1)] }+ = - ela) (z-1)6(2) |? > 0 (2.8) 


where G(z) is the z-transfer function of the linear plant and q is 
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any positive or negative real number. 


Inequality (2.8) is also applicable if the nonlinearity 
Satisties condition, (2./) instead of (2.6) for the case ian. which 


ky = -o and ky = k'+<+o, The constant q in this caseris restuicted 


to be non-negative. However, if ky = -k' and ky = +o, then inequality 


(2.8) is no longer valid and the sufficient condition of system stabil- 


Tiveromiiemsatic lactone OlmtinesciECLe |z| = 1 of the following in- 
equality: 
: eG i ! 2 t 
Re {G(z) [1g] be * <4 (FC (2) ak (2.9) 
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where q is restricted to nonnegative values. 

Jury and Lee showed that Tsypkin's condition (2.2) can be 
obtained as a special case of their results. If no constraint is 
puts,on thewslope of the nonlinearity, i.e. k' = ~, then q must be set 
equal to zero Anveither) (228)s0r (2.9), andi ins either: case (2:2) 
results immediately. 


In two later perverse ‘ 


Jury and Lee extended their results 
to apply to multivariable systems and pulse width modulated eye tere 
Another very interesting work in the study of stability of 
nonlinear discrete-data systems is the work of pelt mehae. He treated 
the problem using an entirely different approach. He investigated the 
system in the time domain. The system is represented by its state 
Space equations. The nonlinear element, which may be time dependent, 
is assumed to satisfy conditions (2.1). By using a similarity trans- 
formation and utilizing the contraction mapping principle, he obtained 
a larger gain sector than that obtained by Tsypkin using his frequency 


domain criterion. This method, like other methods, gives only suf- 


facrent conditions hor “asymptotic stability in-the-larze. 


2.3 Bounded Input-Bounded Output Stability of Nonlinear Discrete- 


Data Systems 
Although it has been proved that for linear systems asymptotic 
ere: ; : Bee ISRAEL 5 
stability implies BIBO stability ,» it has been shown, by counter 
soni that this is not, in general, true for nonlinear systems, 
Furthermore, a nonlinear system can exhibit local BIBO stability even 
though it may not be BIBO stable in a global sense. 


The efforts made in the study of BIBO stability of nonlinear 


control systems were first directed to the case of continuous-data 
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systems. Sandbere and Bergen et epee using two different approaches, 

showed that the satisfaction of the V.M. Popov theorem is also a suf- 

ficient condition so that the. system considered is BIBO stable. 
Nonlinear discrete-data systems were considered by Iwens and 


ells 


Bergen They investigated the same class of systems considered 


by Tsypkin and Jury and Lee. By working in the frequency domain, they 
showed that the satisfaction of inequalities (2.8) and(2. 9) 1s not 
Gniyesurtieient for LConiGagieae stability “in-the-large, but also suf-— 
ficient for establishing BIBO stability for the classes of systems 

they apply to. They also showed, as a special case, that Tsypkin's 
criterion (2.2) also establishes BIBO stability for the class of systems 
it is valid for, even if the nonlinearity is time-varying, provided it 
Seaeioireswconartion: (2.1 D) ror all time. 

An entirely different method, iterative in nature, was 
developed by Suan for investigating the problem of BIBO stability 
of nonlinear continuous control systems. Essentially the system 
response, x(t), is described in terms of the system input, y(t), in 


the form of the Volterra series: 


co 


iH h(t-t)y(t)dt + igi hy(t-t, 5 t-t,)y(t,)y (ty) dt dt, 


CE) = 
SoG 4r fo f h(t-ty> t-Ty see st-t )y(t,)y(t,).-.y(t,)dt dt. 
—0O 7 —0o 
sates C2nL0) 
: th P : 
where h (t-T), t-Ty5+++,t-t) is the n order kernel and is given by: 
i 
h (t-t)> t-Ty 9s -,t-t) = [ h(t-t)h(t-t,)..-h(t-t)dt ee arate WP) 


—_-o 


dt 


~~ — 
a! a aa 
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Tees nit eis the impulse Gh of the linear part of the 
system. The region of convergence: of the series. (2.10) Bee 
the region of BIBO stability of the system under consideration. 

ahtekeeee and Christensen and Tene’ obtained the same 
Series solution (2.10) using functional analysis. The series is 
generated by a mapping which maps the space of continuous functions 
potOmtroe lf: the input and outputsor the system beine selements, of 
that space. By use of the Banach fixed point eieomenc nies the series 
is shown to be convergent, and consequently the system is BIBO stable, 
. in the region where the mapping is a contraction. The advantage of 
this method is that the series solution generated by the mapping 
representing the system is unique in the region in which the mapping 
ds) ay contraction. 

In their work Christensen and Trorree included the effect 
of non-zero initial conditions in the series solution. By setting 
the input to zero they showed that the een where fee mapping 
representing the system in this case is a contraction is also the 
region in which the system is asymptotically stable. 

As far as nonlinear discrete-data systems are concerned, the 
first attempt made to use this approach for the study of these systems 
is due to Rivera. although he did not consider the question of con- 
vergence of the resulting series. He represented the response of the 


System in terms of its input in the form of a discrete; Volterra series 


analogous to (2.10): 
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where x(m) is the output of the system at t=mT, T being the sampling 
period, and hy Ccyo++s sk) is the nth order kernel. Alper used multi- 
dimensional Z-transform to determine these higher order kernels. 

The most extensive work done so far in the field of BIBO 


= He 


stability of nonlinear discrete-data systems is due to Rao 
used a functional representation for the system and by utilizing the 
Banach fixed point theorem he obtained the region of BIBO stability. 
This region is the region in which the Banach fixed point theorem is 
satisfied. Two classes of systems were considered. The first con- 
sists of a nonlinear element located in the feedback path of an 
otherwise linear system. The nonlinearity is assumed to be described 


by: 


b(x) = kx + is ec (2.13) 


where ki k, are constants, and N a finite integer. In the second 
class the nonlinearity is located in the error path rather than the 
feedback path. 

Utilizing the same principle Rao considered BIBO stability 
of systems with slope restricted nonlinearities. Although the results 
he obtained in this case are quite conservative compared to those 
obtained by Jury and Lee, which also apply to BIBO Lyeahpulbiged WS. his 
method has the advantage that the actual solution of the system can be 
generated by a process of successive approximations, starting from the 
solution of the linear system as a first approximation and solving the 
system equation by iteration. The convergence of the solution is 


F ; : F peel 24 
assured in the region in which the mapping is a contraction ra es 
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CHAPTER IIL 


ANALYSIS OF LINEAR DISCRETE-DATA SYSTEMS 


WITH NON-ZERO INITIAL CONDITIONS 


Sole introduction 

In this chapter linear discrete-data systems will be con- 
sidered. Two types of systems will be analysed. The first type has 
no hold olepeitii while in the second system a zero order hold device 
psy included.) Neu—zero initial’ conditions aresincluded in a seneral 
expression forthe solution of the system. BIBO stability as well 


as asymptotic stability will also be discussed. 


3.2 Difference Equation Representation of Linear Discrete-Data 


Systems” *, 

While ordinary differential equations, or vector differential 
equations, successfully describe the dynamic behaviour of continuous- 
data systems, they can not describe the behaviour of discrete-data 
systems since they encounter derivatives of signals occuring within a 
given system. These derivatives are not, in general, defined for a ‘ 
pulse signal in a discrete-data system. 

However, discrete-data systems can be conveniently represent- 
ed by ordinary difference equations, or equivalently by vector difference 
equations. These equations relate the values of system variables at 
different sampling instants. 

Discrete-data systems can be considered to be derived from 


continuous-data systems by considering the variables of such systems 


only at discrete instants of time. This can be achieved by the sampling 
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operation, Fig. 3.1. Most discrete-data systems belong to this type 


of systems. 
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Fig. SL. (b) 


Fig. 3.1 Development of discrete-data 
systems from continuous-data systems 
(a) original continuous system 

(b) discrete-data version 


To show how to obtain the difference equation representing 
the discrete-data system from its continuous-data version, it is most 
convenient to consider a second order system as an illustrative example. 
However, the method is applicable to systems of arbitrary order. 

Consider the continuous time-invariant driven system des- 
cribed by the differential equation: 

x(t) + a X(t) PAE Ge) oennce) (3.1) 
where x(t), x(t) and x(t) is the output of the system and its derivatives, 


¥Ctjeis the input, and ays by» and c, are constants. 
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The corresponding discrete-data system is obtained by con- 
Sidering the quantities in (3.1) at t=kT, T is the sampling period, 
and k=0, 1, 2,...,°. Furthermore, the derivatives can be approximat- 


ed by the following expressions bearing in mind that T is usually 


small: 
Lee ALLE x Co) (3.2) 
Rew a - x(t) 
By making use of (3.2) one gets: 
Beh x(t+2T) - oe + x(t) (3.3) 
ce 


BysSubStLtutions in.(3.1) fwemtiind: 


x[ Cera) T)-2x[ (eth) T)rxlet] 4, xh Cette) oy sper] = c.y[kT] 
1 ih 


mG 1 Ai 
or 
xii) Tleteax| (chi) Tl] bx[kT]) = cy[ kl] (324) 
where a= aT wis 
b= cae = a,T ap 
a= ct 


Equation (3.4) is the difference equation representing the discrete 
VerStOn. OLesystem (3.1). 
Now we consider an nth order linear discrete-data system 
which can be described by an nth order ordinary difference equation: 
aux (itn e+ a oe Ccin—=1 cP bea x (KP) ax (i ay (i) G35) 
n n-1 Ue O 
where a. isea constant, 1=0, 13 25.005) ns. (Kine1) ee tandsetorex( (ktn—i)g)|, 


the value of the output at ta(ktn-i)T and k=0,1,2 ...,%™. 
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System (3.5) can be represented by a block diagram involving 
time delay units, adders and amplifiers. This is shown in Fig. Boe es 
Tesi smclcarstromeiie mma. cethat thesettect of@the inputiwilt 
appear at the output only after n sampling periods. In other words, 
the impulse response of system (3.5) is equal to zero for the first 
n sampling instants. This result can be checked by taking the Z-trans- 


form of equation (3.5). Assuming zero initial conditions one gets: 


ence X(z) + Pe ae X@) eee. ok a,z CRY =e ay X(z) = Y(z) 


which gives 


1) 6 eS ee 


2) @ sp ~l Z, qP 5n0 TF GL Bo a 
n n-L il O 


where X(z) and Y(z) are the Z-transforms of x(k) and y(k) respectively. 


The system z-transfer function is thus given by: 


G(z) ey SA as ae ee ar errs (30) 
ae tral Z a ono oF 2 Bw @ 
n n-L iL O 


The discrete impulse response can be obtained by inverse 
Z-transforming (3.6). Using the power series method it is evident 
Etatetne va lucsoLethewinitpampulse response at t-O0 jn.) 215... (n-1)T 
Sie 

(0) = g(T) = wees. .. = g((n-1)T) = 0 (3.7) 
which confirms the result arrived at from block diagram considerations. 

At this pointeie is very important to note that tThestype or 
systems considered in this section has no hold device. Systems with 


zero order hold devices will be considered at the end of this chapter. 
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3.3 Response of Linear Discrete-Data Systems With Non-Zero 
ead a Ae S2 
Initial Conditions 
In this section we consider a linear discrete-data system 
with no hold device. The system is assumed to be time-invariant and 
of arbitrary order. The difference equation representing the system 
Tseeiven by equation (355) .)vize: 


a_x(ktn) +- ays (etn-1) EP Gas Se a,x (kt1) + a x (x) = y(k) (358) 


where x(ktm-i) ig the value of the system response at t=(k+n-i)T;T is the 
sampling period; y(k) is a bounded driving sampled function; as sea 
COnsStanleicO..),.e6 ,8n neice the order or the system, and k=0, 1, 2, 
++.,”, representing the sampling instants. It is assumed that y(k) is 
ZeromL_oenerativerk® 


By taking the Z-transform of equation (3.8) one gets: 


a 2™[X(z) - x(0) - x(1)z >... -x(n-1) 27] 
+a 2 7[X(z) - x(0) - x(L)2 >... -x(n-2)27*] 
at are 
+ a,2[X(z) - x(0)] + aX (2) y NEGay) 
Rearranging we find 1 
[eee ar feaa ae oA a4z + a. X(z) 
= ¥(2) + [az +a, 2 > 4... + ayzl x(0) 
ae het zo - + a Beg aPogad a,z] ail) 
meres pe ee ws [a2 + a__42! x(n-2) 
+ [a2] x(n-1) (onD 
whererx (0). x(L)5 ue Nob sare the sinitialecondi tions: 


Now we define the initial vector x(o) as: 


I 


Pia) aa ba, a Ge 
iE 


(0) 


2.0) 


[x, (0) x, (0) Ao roe eee x (0) 
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Also let L(z) define the delay operator multiplying X(z) in (3.9). It 


is given by: 


n 
CCA te ge in ee Zz +, wee aos + a 
n 1 fe) 


Making suse Of (3,10) and (3.11) equation (3.9) becomes: 


Ut n-1 n-2 
X(z) ==> 
(z) mS [z(a_z + a2 ie eee a,) x, (0) 
n-2 n=3 
4. 
z(a_z + a2. so Per ea ay) x, (0) 
so, Me ee z(a_z + a__) x16 


4 z(a_) x (0) aC) | 


Coit) 


3 12) 


From (3.12) it can be seen that the coefficient multiplying 


x, (0), Lowel. +) etl canepDe: WrLtren was * 


i n-i n-i-1 
A ies zlaz + a2 ee at a] 





DUbSEbeubingum( al S)minton.(SalAmone gets: 


Nn 
et 3 B, (2) x, (0) + ey Y(z) 
=a 
* | 
= Yi B,(z) x,(0) + G(z) Y(2) 
Es 


G(z) = a defines the z-transfer function of system Gaon 


Now inverse z-transforming (3.14) we get: 


n k 
x(k) =)" by (e) x,€0) + 8G) yr) 
i=l 720 
where 
g(4) = 2 71G(2)] = ge DP Gz) 2 ae 
"a 


and 


ih v, k-1 
b, (k) =e 7, [B, (z)] = ona fp B, (2) 4 dz 





* See the appendix. 
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(3.15) 
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Defining <u, v > as the inner product of the two vectors 


u and v, equation (3.15) can be written as: 


k 
x(k) = <b(k), x(0) > + Sg) yGe3) (3.18) 
5=0 
where 
be at 
b(k) [b, (k) by (k) Pan Aya Aa ts bo] (3.19) 
Equation (3.18) gives the response of the system (3.8) at 
Lae, EO a ewe 


The second term in (3.18) is the convolution summation 


between the sampled input, y*(t), and the discrete impulse response 


of the system, g*(t). Here we have: 


k k 
Dy g(j) y(k-j) = y g(k-j) y(4) 
4=0 ae 
SPS eLeieiie = (3.20) eIntoeoo. | 6) ewe: eet: 
k 
x(k) = < b(k), x(0) > + by g(k-j) y(j) 
He0 


The discrete response of the system is then given by: 


co k 
x(t) = » [<b(k), x(0)> + (4) y(k-§)] 6(t-kT) 
k=0 j=0 
2 k 
= [<b(k), x(0)> + yy g(k-j) y(4)] 6(t-kT) 
k=0 4=0 
where 6(t-kT) is given by: 
= ‘ih for t=kT 


6 (t-kT) 


= 0 for tekr 


(32720) 


Cea) 


(32222) 


(3.226) 


(sms) 
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3.4. Ihe Functional Representation’>’~". 


The dynamic behaviour of any system, linear or nonlinear, 
is usually described by a differential or integral equation, or their 
discrete version if, the fs Rem considered is a discrete-data system. 
However, only in the case of linear systems can one, in general, 
solve the system equation explicitly with relative ease. In contrast, 
in the case of nonlinear systems, it is very difficult, and usually 
impossible, to solve the system equations explicitly. For this 
reason it is very advantageous to use functional analysis. 

1TRLS approachwutilizes as its tools the concept of abstract 
Spaces and operations on such spaces. The input and output a ieneneon 
of the system are identified as elements in these spaces with the 
system considered as a mapping, or a transformation, of the input 
space into the output space. 

An abstract Banach space is suitable for studying discrete- 
data systems. One uses such a space when applying the contraction 
‘mapping theorem to system analysis. The following definitions will 
clarify what is meant by a Banach space. 

Definition 1: A collection of elements together with a 
certain structure of relations between elements or of rules of manipula- 
tion and combination, the whole supporting a mathematical development, 
is often called a space. 

DeSindtion 2: A space, % is Said to be a) linear, wvectonu, 
space if addition and scalar multiplication are defined on X satisfy- 
ing the commutative, associative, and distributive laws. The scalar 
multiplication is related to some associated field of scalars, the 


real or the complex field, usually denoted by 7. 
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Definition 32) Afspace Xs said to be ‘a metric space, 1£ 


with every pair of points x of X is associated a distance function 


es) 


or metric EX, Xy) SuUChme are 


(a) Ax, » x5) S15 TEP f X, 

(b) d(x), x,) = 0, if and only if x, = x, 

(c) d(x), x) = d(x, X,) 

(a) dx, ONG alten Pay ce Witenes ee y) eee ora eee 


3 
Definition 4: A norm on a linear space X is a real valued 
function, whose value at x is denoted by ||x||, with the properties: 
(a) [|x|] > 0 
Gym lac ZO it x 80 
Co) [xy txt] < Ugll + HH 
Gy esp: ey) AMES ssepaewkl 2a ¢, 
The space X with this norm is called a normed linear space. 


[o-e) 
Definition 5: A sequence ie in a metric space is called 


0 
a Cauchy sequence if Ax, x) > QO-as mi andn > =, 

Definition 6: A metric space X is said to be complete if 
every Cauchy sequence of elements of X has a limit in X. 

Definition 7: A Banach space is a normed linear space which 


is also a complete metric space with respect to the metric Ax, » Xx) = 


|x, ae induced by the norm. 


all 
Definition 8: Let X and Y be two linear spaces with the same 


field. of scalars, ands letyA.e we Y be autunction detinedponsxewacn 
Values in Y. A 16 called a linear operaton from xX into; wait ae has the 
following properties: 

(a) ACs, + x5) = A(x,) + A(x,) for any x9 x5 eX. 


(b) A(ax) = aA(x) for any xeX and any scalar a. 
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Now to obtain the functional representation of system (3.8) 
we begin with either equation (3.22a) or (3.22b). Expanding equation 


(3.22a) we get: 


h(t) = 0 < p(k), x(0) > S(t - kT) 
K=0 ap 
. (J) y(k-3)] 6(t-kT) 
k=0  4=0 
= » < b(k) 8(t - kT), x(0) > 
k=0 foe) k 
+ Ys (4) yCk-j)] 6(t-kr) (3.24) 
k=0 5=0 
Now 
fee) (ee) nN 
y < b(k) 6(t-kT), x(0) >= » YP, oCE-kT) x, (9) 
k=0 k=0 i=l 
n (oe) 
C ¥ x, (0) y b, (k) 6 (t-kT) 
i=l k=0 
n 
= Ni x, (0) DaiGc) =e (tne x(0) er 
=i 


Substututine tise back into (3-24) "one “gets: 


co k 
xe(t) = < d¥(E), (0) > + YE ) aGi) yCkrd)] 6(t-eT) 
k=0 j=0 


= £*(t) + B*y*(t) (Gr25) 
where £*(t) = < b¥*(t), x(0) > (226) 


and ®* is an operator defined by 


ee k 
gey*(t) = ‘a is (4) yCe-4)) 6k) (3.27a) 
k=0 4j=0 


Of, Starting with equation (3422p), 
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0 k 
Bxy*(t) = » » g(k-j) y(j)] 6(t-kT) (3.276) 
k=0 4j=0 


It is easy to prove that 8* is a linear operator. 

Now the input y*(t) and the output x*(t) are ro valued 
functions of time defined only atethe discrete instants t = kT, 
Keen oes se, ©. lhey earesthus sequences in a sequence space, 
say S, which is a linear space. Addition and scalar multiplication 
in this space are as follows: 


Addition: If x* = io he oi 
LteLOn cx aaa and y iva! oetnen 


zk = x* + y*, where z* ={z }~ 
2 n n=0’ 


Multiplication: ax* = a{x }” | = {ax }- , ac#. Here F 
n n= n’ n=0 


OA ES Ge 5 
n n Tn 


tom themiaeldsorereal scalars, 

If furthermore, y*(t) and x*(t) are bounded, they belong 
to the space of bounded sequences, say F, which is a subspace of S. 
F is a Banach space if we define the norm: 


Pix*®() | [= sup” “[x*(#),] 1p xxer. (3.28) 


—0< f{ <oo 
Now equation (3.25) can be written as 
x*(t) = A® y*(t) (3.29) 
where A* is an operator defined by (3.25). A* is a mapping of S into 
itself in the general case, while for the special case where y*(t) and 
x*(t) are bounded, A* maps F into F. 
Sey oteapility of linear Discrete-Data Systems>>. 

In this section asymptotic stability as well as BIBO stability 
of the linear discrete-data system (3.8) will be described. As a pre- 
liminary development we will discuss further the coefficients B. (2), 
and consequently b, (k), inal en eGekined by eduatdonsa(3.13) 


and (3.17) respectively. 
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B, (2) can be rewritten as: 


tL n-itl n-i 
Loy [az + a_12 ce eae a2] 


Il 


B. (2) 


1 Be n=-1 
G(z)[a.z inoue; Z SCRE. as a, 2] (3.30) 


Now, according to (3.7) we have: 


m-1 
& eg! 4 
Pee Caine! | 25 Gamez Sy 2 
J=0 
ae 2eG(Z). ml = lee y, ae Si (8.031) 
or 
Se oe i) m 
eeCerml je 2 “(2 GG) mw, Uae nt (332) 
Inverse: z—transiorming ™(3,.30) and ubiligime (3.31) and (3732) 
we find: 
be (Ct) = a s*(tt(n-it])T] + a e*[tt(n-i)T] 
aP 264 SF a,8*(ttT) (35.35) 
or 
b, Ck) = a e(ktn-itl) + a _ 18 (ktn-i) + ae oe ag (k+1) (3534) 


Se me wee ee, oe wis la tye lbs 


Sips seu e nee. Oo ine Ol G3..l5) sweet ind: 


n 
a(k) o> \ [a e(ktn-itl) a a _18(ktn-i) age BP a, g(k+1) |] x, (0) 
i=1 
a ye g(j) y(k-3) (3°35) 
Few 


Se saleeBLBOnS lability. 


By taking the absolute value of both sides of equation (3.35) 


OnemaeteSes 
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Ix(k)| < ») | ag (ktn-it1) PEW yelC gaint) © eons a,g(kt1)||x, (0) | 
ee 
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le(3)||y(e-3) | 


(aay 
w [~]= 
° 


n 
<) max Ja,| |e(ein-it1) + g(ketn-i) +... + gC] |x, (0) | 


Cc 


+ sup |y(k)| y les) | 


0 sis 


j=0 
= 7 re aes [1£(-it1) sup |g(5)|]1]x, 0d] + |]y*¢e) || sy |e(3)| 
G1 051,52 _ 
3=0 
<) max |a,|. max |x, (0) easupee (4) e né 
l<j<n l<i<n O<4s2 
Ge | » le(3)| (3. 36) 
j=0 
Since the input is bounded and |]y*(t)|| is finite, the 


ReivoemrOrminequalbity (3.06) is bounded it ; 


oO 


Ye leG)| <@ £43537) 


ENUSmownicem(Gas0) Usetrue Tor all kk. (5.3/7) 18 a sure 


ficient condition that ||x*(t)|| is bounded. Also it has been 
proved that (3.37) is also a necessary condition for BIBO aaaper aig 
Notes that’ (3937) *inmplies=thats 
lin |e(5)|e= 0 (Gy 3D) 
jr 
Dae eens yMplLOtuc oe tabi lity 
For the purpose of studying asymptotic stability, we let 
the inputsok the system y= (2), ,equall toezero ing (3.35) stoeget. 
n 
x(k) = 5 [a g(ktn-itl) + a _18(ktn-t) “2 Geb as a, 8(k+1) ]x, (0) 
i=1 
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Now asymptotic stability requires: 


Vim |x(ky| =_0 (3.40) 


ke 


From (3.39) we finda. 


nN 





PGs ‘3 {Ja | |e(ktn-it1) | “| [a__; | le(etn-4) | +...+/a,||g(kt1)[}|x, (0) | 
i=1 
(a1. 
Clearly, since a. and x, (0) LG oy 1h care finitesand sare not 


allezero,~nence, 11 (3.37) holds ,tnen so does .(3<38) and the R:H.S. of 
(eal jewil Metend to zerojassk tends to infinity. Therefore (3.3/7) is 
also a necessary and sufficient condition for asymptotic stability of 
linear discrete-data systems. 

3.6 Systems With Zero Order Hold. 

For the purpose of investigating linear discrete-data systems 
which include zero order hold (ZOH) devices, it is again preferable to 
represent the system by a differenceequation. However, in this case 
the approach of section 3.2 can not be used since the differential 
equation representing the plant together with the ZOH is rather com- 
plicated. lt seems logical™ that’ to:determine the required) difference 
equation, one can make use of the z-transfer function of the whole 
evetenin aoe 

To illustrate this we consider a system whose overall z- 
transfer function together with the ZOH is given as a ratio of two 
polynomials in z. The number of its zeros is at least one less than 
the number of its poles. This is usually the case in most physical 


systems. The impulse response of such a system will be zero for t=0. 


Imeother words 2@(kl) = s0etonsk=0. 
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Mathematically: 
m m-—-1 
+c Zz a2 bop ae © 
X(z) <n m-L fe) 
G — = 
€2) YeCz) n n-l erie 
ay we abe) Z aera Hur ees) 
n n-L O 


where X(z) and Y(z) are the z-transforms of the output and input 
respectively, and m < n. 

sinee the zg-transier function G(z) is’ independent of initial 
conditions, cross multiplying (3.42), inverse z-transforming, and 
SAlcineeinto sonsideration relation? (3.32), an gets: 


ax (kn) + ax (ketn-1) ce are a x(k) 
= cy (+m) =e c py (etal) eee eke coy (k) (3743) 


Equation (3.43) is the difference equation representing the system. 
Now to obtain the response of the system when initial condi- 


Evens are not zero, we take the z-transform of equation (3.43). 


Hence on. 


n m 


ioe 2 Cede 7 B, (2) x, (0) - ». F, (2) y, (0) (3.44) 


i=l i=l 


where the initial vectors x(0) and y(0) are given by: 


OEE (ORCL) ceearex(n-1) la 
> EG) GQ) os 2 5 Or (3.45) 
and 
y(0) = [y(0) y) . . . . y(n) 1" 
S(O) et) 5 oo 5 Or (3.46) 


B (2) is (elven bys Cowboy), and F, (2) is given by: 


=: ey 
F(z) -= S ey eitee Se aie nce C5] (3247) 
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By following the same procedure as in section 3.3 we get: 


k 
x(k) = < b(K), 2(0)> - <£(k), y(0) > +) aCsdy(k-3) (3.48) 
4=0 


moeresp(k).) (kc) andme Giomaremcetinedibystelations similar to (3.17), 
(ano eand .G3.16).. 


The discrete response of system (3.43) is then given by: 


oo k 
xe(t) =) [eb(ke) ,(0)> - <£(k) ,y(0)> + Y° gC5) yCk-J) 15(e-k) 
io | ie (3.49a) 
k 
= [<b(k) ,x(0)> - <£(k) ,y(0)> + » g (k-j)y (4) 16 (t-kT) 
ta | de (3.49b) 


where 6(t-kT) is the delta function defined by (3.23). 
As far as system stability is concerned, an argument like 
that of section 3.5 can be used to show that conditions (3.37) and 


(Sees) aresalso applicable for systems with, ZOH devices. 
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CHAPTER IV 


ANALYSIS OF NONLINEAR DISCRETE-DATA SYSTEMS 


WITH NON-ZERO INITIAL CONDITIONS 


4,1 Introduction 

In fits: euapter we will consider time invariant, nonlinear 
discrete-data systems. It will be shown that the response of these 
systems can be represented in the form of the discrete Volterra 
series. The series is developed by solving the system equation by 
iteration. The convergence of the series, and consequently system 
Stability, is investigated by the aid of the Banach fixed point 
theorem. 

Systems without hold circuits will be discussed first, and 
then the results will be extended to systems including ZOH devices. 
4,2 Description of the System 

The configuration of the nonlinear discrete-data system 
under investigation is shown in Fig. 4.1. It is a single-loop nega— 
tive feedback system with one nonlinear element N located in the 
feedback path. The forward path consists of a sampler followed by 
a linear plant Ge. The sampling is uniform and of period T. 

x(t), yGt) and e(t) denote the “output; input and error signals) res— 
pectively; they are continuous real valued functions of time. 


x*(t), y*(t) and e*(t) are the corresponding sampled functions. 
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Fig. 4.1 Nonlinear Discrete-Data System 
Under Investigation 


The nonlinear gain function $(x) is assumed to be continuous 


in its argument, memoryless, and further satisfies the condition 


¢(0) = 0. These assumptions allow its representation in the form: 
M 
d(x) = yx + s yx" CAB iL) 
i 
i1=2 
yx + 6,09) (4.2) 


where Yy> Yy are constants, and M a finite integer. 
It is assumed further that the corresponding linearized 
ey ecen obtained by neglecting 1) is BLIBOMstable. 
4.3 Response of Nonlinear Discrete-Data Systems Containing 
No old) Devices 
The system considered in this section has no hold device 
included in the linear plant Ola The forward path can thus be re- 


presented by a difference equation of the form (3.8): 
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ax (k+tn) +a x(ktn-1) Aca oee eet a x(k) = e(k) (73) 


where x(k) and e(k) are the output and the error signals measured at 
t=kT respectively. 
But 

OMB ny (Kh) = OCR (k)) (4.4) 
then, substituting (4o4) i ben (es )e-and making use of Cia2)) cnesgets: 


a_x(ktn) + a yet“) SP pa Be (a +y) x(k) 


= yk) - $, (x(k)) (4.5) 


Equation (4.5) represents the nonlinear discrete-data system under 
investigation. 
Notememat if $ (x(k) is) sebwequal to zero in) @.)), we cet: 


a x (kn) + a px (etn-1) oe Gobi SF (a ty) x(k) = y(k) (4.6) 


whiten is the equation of *the linearized system corresponding to the 


nonlinear system (4.0) © The’ solution of (4°96)01s thus. givem by: 


© k : , 
ve(t) = ye [< b(k), x(0) > + y ai) pO OCS) 
K=0 | 5=0 
or 
ve(t) = £*(t) + B*y*(t) (4.8) 
where | 
PGE)! Fe poled in. SG hye (4.9) 


B* is the convolution summation operator. ‘It is defined by: 


co k 

Bey*(t) = s a eG )< fe) 8c a) (4.10) 
k=0 4=0 

BG 2 Ie) 1 need (4.12) 


4 = 
L(z) suave wa) fe ©. Hen ta, 2) (acty)) (4.12) 
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New wectorsb(k)eisidefined by equations (3.19), (3.17) and (3,13): 
Also the initial vector (0) bo defined. py equation (3,10). 2@(]) is 
the impulse response of the linearized system evaluated at t=jT. 
Now, we consider the nonlinear system (4.5), and we let 
e,(k) = y(k) - 9, (x(k)) » (4.13) 
Then, equation (4.5) becomes: 


a_x(ktn) + ax (ctn~1) +e, PoP (a sty) lc) a= e, (kx) (4714) 


which has the same form as equation (4.6) with e, (x) taking the place 


Obey Cc uinus, (4. 14)has arsolution: 


© k 
xA(t) = Yo [< b(k), x€0)> +) 2C4) ey (k-3)] 6Ct-kT) (4-15) 
k=0 j=0 
or 
seep) ee Sac Ge) ae B¥e,*(t) : (4216) 
Since 8* is a linear operator, expanding (4.16) using (4.13) 
we get: 
ae eet) ery k(t) > BO, (x(k) (4.17) 
= v¥(t) — BOT (x(k)) (4.18) 


where v*(t) is the solution of the corresponding linearized system 
given by equation (4.8). 
Now. we substitute ain (4.17) and (4.18) with) the valucgor 


91) obtained from equations (4.1) and (4.2) to get: 


M 

h(t) = ERC) + BayR(E) - BRL YO vy (DI (4.19) 
* i=2 | 

= v¥(t) - 8*| ). yx (t)] , (4.20) 


i=2 


Withetie inpue y(t) fixed  equarion (419 ore (4.20 ledetimes sal rela— 
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x*(t) = y x*(t) (4, 2) 


where y is a nonlinear operator defined by: 


M 
U x(t) = £¥(t) + BAy#(t) - B*L ae (4.22) 
i=2 
H +k 
2 GVS Ty ug CN (4.23) 
i=2 | 


Equation (4.20) can be solved for x*(t) by a method of 
successive approximations starting with the solution of the linearized 
Sysctemev-(t) as a first approximation. glhis is carried out as follows. 

Without loss of generality, we assume that the nonlinearity 
is given by: 

ly S Ay Be eS (4.24) 
which is obtained from (4.1) by letting y yA, Ya= e>0, and 
y,=0; Loe DR ME), Mtie re 5 aha 

The nonlinear difference equation representing the nonlinear 
system in this case is obtained by substituting (4.24) into (4.5) to 
get: 

ax (ktn) + a _yx(ktn-1) SP eth PR (a +1) x(k) = y(k) - ers (425) 


The solution of the system is, thus: 


3% 


x*(t) = £*(t) + B¥y*(t) -. eB* x” (t) (4.26) 
B yet )e tops LACE) (4.27) 
The operator U is also given by: 
Uxs(t)M=NEs(t) Ses ys(e)e= ese x? (t) (4.28) 
= vA(t) - Bt x” (t) a (4.29) 


To carry out the process of successive approximations, it is 
more convenient to rewrite equation (4.27) for the output at the end 


of the kth sampling period, thus; 
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3 
ple) = =v (ke) ge Bee (1) (4.30) 
where 8 represents the convolution summation operator for that samp1- 


ing instant. 8 is defined by the relation: 
fe ; 
Bu(k) =) g(k-4) ud) (4.31) 
a= 
where g(n) represents the impulse response of the linearized system 
at the end of the nth sampling period. Utilizing (4.31) in (4.30) 


Onew2e es: 
k 


E(hyee (ey ay ee Ck) (4.32) 
5=0 


Now, the solution v(k) of the linearized system will be the 
first approximation xo (k) tae 
x Co) = v(k) (25 33) 
Subse mcutine (4,33) into theek.H.og,o0f. (4.32) we obtain the second 


approxima tLonm LOL, x(k), Viz. x, (x) : 
k 
x, (k) = v(k) - ¢ ye E(t) (4.34) 
5-0 


Burther isubstitution of (4234) into the R.H.S. of (4.32) results in 


the third approximation: 


k j 
x(k) = vl) = ©) gg) (VG) - © Y > gG-d) WGI 4.35) 
4j=0 jae 


The second term in (4.35) can be expanded in the following way: 
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j 5 
2 
225 SOM) ap Gir 
j,-9 
3 3 3 3 
a Il AC ak ae Ge Py) 
j,=0 
k k j 
=) a(k-§) via) - 3¢ o 2 a(k-J)e(-5,) VG) VG) 
au) j=0 j,=0 
k j j 
SL), Aeaneeaneey von ee) clan 
aU! aly 
k J j j 
ee bs bs 2, [g(k-j) 8(4-3,)8G-5,)8G-3 ae (3 pe Gi5)" 
Pj 31 (4. 36) 
But since as n) = 0 for all positive integers n, we can write: 
i ec NeG-3,) VG) VG) 


eu) tee 


k 
ue g(k-$)8(3-4,) vw G) WG) (4.37) 


iM 


Similarly for the other terms. 


Substituting (4.36) into (4.35) making use of (4.3/7) we obtain: 


k k k 
x(k) = vo) -2 S7 ates) VG) #327) DD ate DeG-apw GWG) 
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k k 
-36) ees BOD 8-58.) VO" GPW Gy) 
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as s 3 3) 3 
te », y ) | a(k-s)a(5-4)8G-4,)8G-Jg)¥ Gv Gy)¥ Ga) 


(4,38) 
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To obtain higher order approximations we must then sub- 


etitute equation (4.38) into the R.H.S. of (4.32). Continuing this 
process, we obtain the solution for the nonlinear discrete-data 


system (4.25) in the form of the infinite series: 
k 


x(k) = v(k)i- e » g(k-j) v°(j) 


k 
v 
+3) Yates pw vr G,) 


k k k 
y net. Chiee at eee RO N. 
3¢ », oD g(k-4)8(i-5,)8(i-5,) Vv Gv Gi) 
4 = 


j= 
k k k k 
A 
teh YY YY ete 06-5 e6-3,)8G-1,)V G Pv G,WG,) 


4=0 j4=0 a) syed j,70 


ole 


Sle Gas ere) ot (4.39) 


Note here that if the system is subject to zero initial 
conditions, then f(k) = 0, and the solution of the linearized system 


will be given by: 
k 


ome ys ea GD (4.40) 
j=0 


Substituting (4.40) into (4.39), the series solution of the nonlinear 


Z. 
system will then reduce to the form obtained by Rao 7g 


Gk) = 2 AVP AGN) vc y ye ssyGpyG, )yG,) 


Se ey Pawee oe (4. Ads) 
where g. is defined by: 
B 
k 
: ma “ea i—4 ea Ure 
a=) a(k-j)g(J-5,)e(I-4,)85-4y) (4.42) 


5-0 
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Now by taking the absolute value of both sides of (4.39), 


we get: 


k 
xe] < Iwo] + fel Jeaeay | vcs)? 
j=0 
ora 
2 
+ 3le| ey 3 JeC-slgG-s ply I71vG 1? 
- 0 Ki = 
+ 3/e| ay a (e-§) | 8-5.) | [8G-3.) | [VD [ |v) Plvas1? 
dOn ihe elon? 
Bur age ie Bato 
k 
< sup |v(k)| + Je|[ sup |v(k)]] > Y lete-5)| 
-O<k<e O0<k<o 
a {20 
2 
Tess | Wesup ivik) Ceaa) Pete ) 
[el "{ sup [ve |] ‘y > [ac-3) | !e3-4,)1 
Fe 1 OTin =? 
k k k 


3 di 
ay [ s (k) |] Ce) ee Gal te ade) 
isla sragle SS YY lets) [1eG-5) 15-4, 


A PEE aL ew 
een ae (4043) 


Since (4.43) is true for all KS it canebeRwritten in the form: 





eee] s [Ivace | + Tel Leg] Pivece [1° + sel Hlesll Peco? 
+ 3fel*IIg,I] [vec ]|7 +... 4.44) 

or be 

Hcy] < Tivecoll + ¥ dat Hell [vece i? (4.45) 

n=1 7 

where in (4.44) and (4.45). we have: 

[}v*(t)|| = sup |[v(k)| (4.46) 

O<k<e 


A similar expression holds for | |x*(t) ||. 
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Definition (4.46) is based on the assumption that 
v*(t) = 0 for tae () (4.47) 
Now, if the nonlinear discrete-data system (4.25) is BIBO 
stable®, then ithe R.H.S...0f .(4.45) is convergent. If also the system 
possesses a unique solution”, then (4.45) shows that the system is 


analytic and equation (4.39) is a Volterra sor ea 


4.4 Nonlinear Discrete-Data Systems.with ZOH 

The system to be considered in this section has the same 
configuration of Fig. 4.1. However, in the present case the linear 
plant G, contains a zero order hold device. Thus, the forward path, 
being linear, can be represented by a difference equation of ae 
form (3.43): 


a_x (cin) a a _ yx (crn-1) SM angie a x(k) 
= ce (tm) + ce (etm 1) Se 6G Sh ce) (4.48) 


where m <n for a physical plant; and x(k) and e(k) are the output 
and error signals at t=kT respectively. 
The error signal e(k) is given by: 
eG) 2 snl) ered) | (4.49) 
Substituting (4.49) into (4.48) and making use of (4.2) we find: 


a x(ctn) + a _ yx (ktn-1) am duro ae a pps (ktm) + (a to )x (kt) 


+ (a to) x Cctm-1) Ore ee ear (a +05) x(k) 


m-1L 


= cae Cet) + Cp ey (et D “he aS os ce, 09 (4550) 


where a, = Y4Cq> Lope ke eee C455 1) 


EO ea AGRE om ERE TEND, (4.52) 


i 


Sea section 4.5. 
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Equation (4.50) describes the nonlinear discrete-data system under 


investigation. 


It can be readily seen that equation (4.50) has the form 


C7 o)= thus-its solutton is: 
co k 
x*(t) = y {<b(k) ,x(0)> - <£(k) ,e, (0)> a4 », g(k-j)e, (5) }6(t-kT) 
k=0 4=0 (4053) 


Note here the e, (0) is a known vector, the ith element of which is 


given by 

55, GOV S PAO ans hpi als agi ea (4.54) 
so that 

Br mee OP sh (0))) (4.55) 


Taree een as MaRS MRE (O RO meats CONeTe 
defined by equations (3.45), (3.46), (3.19) respectively. £(k) is 
depinedsby ecelations similarn.to. (3.17) .and,(3419). o,' )) is 
defined by: | 

$5 (8"(0)) = [6,0%,(0)) 6, 6,(0)) -  . 9,01" (4.56) 
and 
1) = 240) © 2), CONIA, x (01 (4.57) 


The discrete impulse response is obtained from: 





g(j) = Z (G(2)) (4.58) 
where 
@ zg ar @ aoe of een te 2 ete ce 
G(z) = a z a a 1 ake ee (4.59) 
By be ae a1 DU NORON on (a ta )z + -.t(a sta) 


Substituting equations (4.52) cand, (4.55)) into, (4.53) we 
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co 


A(t) =) L<b(K) x(0)> - <£00 .y(0)> ~ <£00, 4, x" (0))> 
k=0 
k k 


+) g(k-ayyG) - Y, BKI)44 OCG)) }5(t-kT) (4.60) 
j=0 j=0 


It can be seen that if we neglect the nonlinear element 
o, G9) from (4.60), we obtain the solution, v*(t), of the linearized 


System=etven Dy (3.49) ,~so that we tet: 


© k 
we(t) = ve(t) =) {<£(e) 49, 60" (0))> + D849, GI) 16 Ce-kT) 
ae 150 (4.61) 
or, using functional representation 
x(t) = vA(t) = wk(t) - B*O*(x(k)) (4.62) 
where © 
w(t) = ». <£ (Kk) 16 (x'(0))> 8 (t-kT) (4.63) 


k=0 


and 8* is the convolution summation operator defined by equation (3.27), 


WellZecn bs : 
erur(t) = ) TY g(e-jduCi)} (t-kn) (aes) 
k=0  4=0 
00 k 
oe i ya BC) alka). (ic ldh) Chay) 
k=0 qc 
If, we nowssubstitute; for >, GCc)) into (4.62) one gets: 
M ok 
h(t) = v(t) - we(t) = BRL Y yx (8) ] (4.65) 
{=2 


With the input y*(t) fixed equation (4.65) defines the 
operator equation: 
x*(t) = U x*(t) 


where Uis a nonlinear operator given by: 
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M 

Ux* = vA(t) - w(t) - BAT ¥ ygxt (t)] (4.66) 
i=2 

ll pecdtatLons 


Now if we let Y3 =o, 6-0 seand Mee OSRr es 2 5 echoes 


(4.65) and (4.66) become: 


n 


x(t) = ve(t) - wk(t) - e@x? (t) (4.67) 
and 
; : op 
{| dees = v*(t) - w*(t) =- ERX Ce) C4, 68) 
Equation (4.67) can be solved by iteration, starting) with q*(t) as 
a first approximation, in the manner described in section 4.3 to get 


thes Volterra series: 


k 
x(k) = q(k) - & 3 SOapenriep 
j=0 
; k k 
2 
$37 VY gtede-5) 7G) PG) 
j=0 4,0 
Sn a ee ee (4.69) 
where q(k) = v(k) - w(k) (45/0) 


Note q(k) is a known sequence since both v(k) and w(k) are known. 
Nocewalso nate (4.09) is similar to (4.39) swith q(k)) taking the place 
of v(k), and thus it can be justified in the same way. 
4.5 Convergence of the Discrete Volterra Series and 
Stability of Nonlinear Discrete-Data Systems. 

| In this section we will investigate the stability of the 
nonlinear discrete-data systems considered in this chapter. This will 
be done by considering the convergence of the Volterra series represent- 
ing system response. Convergence will be established by using the 
contraction mapping theorem. The theorem is stated next. The proof 


of the theorem is given in referencé 25. 
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THEOREM Let F be a Banach space, and let A be a mapping 


OmretntO ttself. Met Ss Dea sphere contained in F of centre u:, 
Oo 


and radius d, such that: 


(i) | |Au, - Au 


JA 


rea as | | (4.71) 
for every Uj> Uy & 5S and some A: O<i<1 
(ii) | ]Au, eared lees: (hs ah va (4.72) 
then the mapping A possesses a unique fixed point u in S, that is a 
point for which 
Au = u CARED) 
Condition (i) will be called the contraction condition, and 
condition (ii) athe fixed point condition. 
Now, to illustrate the application of this theorem to the 
investigation of system stability, it is most convenient to consider 
a specific example. 


Illustrative Example 


Consider the nonlinear discrete-data system shown in Fig. 4.2. 


x(t) 





Fig. 4.2 Nonlinear Discrete-Data System 
For Illustrative Example 


j | rl > sth = AT | (i) . 
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The system contains a nonlinear element located in the feedback path 


which is described by the gain function: 


d(x) =xte oe (4.74) 


The forward path of the system consists of a sampler with period 
T=1 second, followed by a zero order hold and a first-order linear 
plant whose transfer function, Cis), is given by: 


ih " 
CS Borers C75) 





INeizatboane pomirinction Of the forward pati is, thus. elven, by: 








Go (2) = Z[ 





| 40.632 
~ 3-0. 368 eID, 
The linearized system is obtained by letting e=0 in (4.74). The over- 


all z-transfer function of the linearized system is given by: 
G (4) 


G(z) = LG, (2) 


iy le Ga) 


clearly, the linearized system is stable since its only pole lies in- 
Sidesthesinitecircle | z| uel, Tel eres 
The difference equation representing the nonlinear system is readily 
obtained as: 

meting 0h 2 G4ex (OI 0463286 oC) = rool (4.78) 
where y(k) is the value of the input function, y(t), at time t=kT=k. 
We note here that the difference equation representing the nonlinear 


system, which contains a ZOH, is of the form (4925) erpathersthan (4.50). 
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This is because the linear plant G(s) is a first order plant. The 


solution of the system is, thus: 
k 


09) = 109) os BACs) Xe) (4.79) 
j=0 


where v(k) is the solution of the linearized system. It is given 
Dy: = 


Tn eo (0) ». ee) (4.80) 
j=0 


where g(k), the impulse response of the linearized system, is given 


by: 
: k-1 
g(k) = ye Residues ofeG(z) z at the 
poles of G(z) Be (4.81) 
0 for k = 0 
i (4.82) 
0 .632(-0.264) <7 sjeiigesal 
and 
fe Zz 
Be ea ca%en 
k 
= (-0.246) ‘ keen0 (42.83) 
Equation (4.80) can be written as: 
v(t) = £*(t) + BAy*(t) SPie 
where 
ce CG) pb (trex (0) (4.85) 
Taking the norm on both sides of (4.84), we get 
Vs = || ve ce) | | = sup | vk) | 
0 <k<co 
< Fe + He Y* (4.86) 
where 
Fx = 


|Jexce)|| = sup [£(k)| =[x(o) (4.87) 


O0<k< 
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us = |[e4l| =) Jace | 
k=0 
= 3 hmen 70m 6) | 
k=1 
=). 632 y (e264). 
k=1 
= 0.86 (4.89) 
Now, equation (4.79) can be written as: 
x*(t) = ve(t) - Pern CL) (4.90) 
=U xe Ct) (4.91) 


To apply She Sanecdiesten condition (4.7/1), we let x# (t) and 
it) be two elements of the sphere ScF, centred at v*(t) and of 
Raamicw a oF (he choice of v*(t) as the ‘centre ‘of S)is justitied on 
the basis that the nonlinear system may be looked upon as a perturbed 
version of the linearized system, so that if it has a fixed point for 
a given y*(t) € F and some initial conditions, then this fixed point 
is likely to be in the neighbourhood 'd' of the solution of the linear- 
ized system for the same input y*(t) e F and same initial conditions. 
Now “the, contraction, condition, is; 


| | Ux; (t) - Ux5(t) || <A | [x4 Ce) - x5 (t) | | Ha << Ve (4.92) 


From equation -(4.90): 


paces ux; (t) || 
= eat? (e) - optx? (t) || 
< el [le*l] [Ce - “1 Ol 
< fel H* [[tcm - FCO] Yo Lge feo? 


ee (4.93) 
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a * * F 
ince x, (Ct), X(t) ess 


ce) LL [keep [Lt lacey || + 


= Ve +d (4.94) 
then: 
ig *(t) be, * pape Z * * 
x] U mE) | < 3 Je] H*(v* + d) | 1x5 Ct) - x7 (t) | | (4.95) 
since | [x5 (e) - my ct) |] = | [x7 Ce) - x5 (t) | » (4.95) is equivalent 
e© 
3 le| Hx(ve + dy? < (4.96) 


We now consider the fixed point condition (4.72), 


|] Uvece) - ve(e)|| < C1 - Add (4.97) 


From equations (4.90) and (4.91), it follows: 
U v*(t) == v*(t) _ eBevs(t) (4.98) 
Transposing v*(t) to the L.H.S. and taking the norms on both sides, 


Onem Zeus: 


seer) en eben s(t) 
< fel [eel] []veco 


= [| Hx yx3 (4.99) 


Hence (4.97) is satisfied if: 


|e| H* ye <eGli= X)d (4.100) 


Races considered the same class of systems, but with zero 
initial conditions. He obtained the same inequalities (4.96) and 
(4.100) with H*Y* taking the place of V*. For the particular example 
considered here, he obtained as a bound for the input, such that the 
nonlinear system is BIBO stable, the following: 

Ye <= 0455 Bi (4,101) 
which corresponds to: 


H*Y* < 0.374 aie (4.102) 
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This means that, for the nonlinear system to be BIBO stable, the 
response of the linearized system must be bounded by the quantity 
TnPCNe wR wiso. OluCs, LOZ), gelm ourwcase, this corresponds to: 

We Pee He YX eno 0 374 je|~* (4.103) 


which gives as a bound for the input: 
ib —l, 
Rie [Oeai4)e| 6 =2Re] 


= [OmO74 pn ~|x(0) |] (4.104) 


If the system is undriven, Y* = 0 and we get: 


1 
78 


agen scO <0. 3745) ch (4.105) 
which is the bound on x(0) such that the nonlinear system is 
asymptotically stable. 
Moree = 0.1. conditions) (4.104) and (4.105) become, respect— 

ively: 

Vem e163) [tees x (0) |] (4.106) 

[x(0):/ <1 41 83 (4.107) 
Clearly it we let x(0) = 0 in (4.104) we get 

Vest lees 15 (4.108) 
which is the result obtained by Reeae 


Now we recall that the response of the nonlinear system 
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ie ie le ae . : 
4 
tee 2D. » y ». [e(k-4)@(5-34)8(5-4,)8(5-54)¥° G Dv Gp) 
= 


Ae OA he 
uh 2 37m 
aa) 
<P nin ect oie eae eMey arty. te (4.109) 
By taking the norms on both sides we found: 
a a 30. 
2S Nes y EP Mesa) i (4.110) 


Hence, since (4.96). and (4.100), ane isatisiied, the 2terative 
series solution found from (4.79) is a unique convergent series. Further- 
more, since the solution can be put into the form (4.110), that series 


: ; : , 24 
is the Volterra series solution for the system considered . 
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CHAPTER V 


RESPONSE OF NONLINEAR DISCRETE-DATA SYSTEMS 


BETWEEN SAMPLING INSTANTS 


Sl in eroduction 


Since the response of sampled-data systems is continuous 
in time, it is very important to know how they behave between sampling 
instants. A system may seem to be stable judged by its response at 
sampling instants, though it may exhibit an unstable response between 
sampling instants. 

In this chapter nonlinear discrete-data systems will be 
considered. Relations that determine the intersampling response of 
the two types of systems considered in chapter IV will be derived. 

It will be shown that these responses are obtained in terms of an 
infinite series. Stability of the systems between sampling instants 
is considered at the same time. 

5.2 Deriviation of Response Equations 

The system considered is described in section 4.2 and is 
shown in Fig. 5.1. The linear plant Go will be defined by its trans- 
fer uch Lon G(s). it may or may not intludevarz0h device. The 
nonlinear element N in the feedback path is described by relations 
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Gc) aie eer Naser (5-1) 
i=2 
= 71% + 9, ) (5.2) 
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It is well ener” that the intersampling response of linear 


discrete-data systems can be determined from the relation: 
k 
(ktm) =) g(ictm-j) (4) Gy) 
4j=0 


where x(ktm) is the response of the system at t = (k+tm)T, O<m<1, 


Pets eunessempling period: y(}) is the input at tejT and (ict) is 


given by: 
vil 
g (k+m) aR [GGed (5.4) 
t=(k+m)T 
where Pet ] represents the operation of taking the inverse 


Laplace-transform. 
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Fig. 5.1 Nonlinear Discrete-Data System For 
The Purpose of Determination of 
Intersampling Response. 


Applying relation (5.3) to the forward path of the system 
Ge eigu lo alewe tind 
k 
x(n) =) g, (ketm-j) (3) 
j=0 


(525) 
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whereve(j|)§is the error signal at t=jf. It is given by: 
Bie ty Ge toh) (5.6) 
Sebotieuting (5.6) into (575), and making use of, (5.2), equation 


(5.5) becomes: 
k 


x(etm) =) g (eem-j)[yG) = yx) ~ 0, 600] (5.7) 
j=0 
Note that if >, 6«G)) is zero, the system will reduce to the linear- 


ized system, whose response will be 


k 
viktm) =) gy (ktm-3) Ly) - 7,91 (5.8) 


4j=0 
The intersampling response of the nonlinear discrete-data 
system is now readily obtained by substituting the response of the 
eyvatem at sampling anatants in the R-H.S. of (5./) 
For the particular nonlinearity considered in chapter IV 
which is: 
Rete (Oa eae) (5.9) 


equation (5.7) becomes: 
k 
x(tm) =)” g (etm-3) fy(§) - x) - ex GDI (5.10) 
j=0 


Now, for a system that does not contain a ZOH device, the response 


at sampling instants is given by equation (4.39), viz.: 
k 


x(k) = vUk) = «€ a ep “in cate (5 ele) 
j=0 


Substituting (5<lL)) into (5.10) one gets: 
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k j 
xem) =) a (kim-siy() - vG) +e Y gg-apwGp- 
i=0 i,=0 
j 
2 eat) =e | Cer Ghp ae 
dies (5.12) 


Expanding equation (5.12) we get: 


k 
x(kctm) = » g,(ktm-j) [y(j) - v(4)] 
j=0 
k j 
eRe ocem-piv@d - Yo sG-sipwayl 
j=0 j,=0 
Fee : g mt) gGai pe (wiG p73 ; Yat j,)8G-j,) 
=O | ‘hie 0 j,=0 
V Hes a (i, )] 
Hee So ae nen ik os (5s) 
k al 
= v(kt+m) - ¢« iy g, (etm-§) Lv (4) ~ iy gG-i)v G1] 
j=0 any 
k 
fa go ys 9 (ktm- ay D eel or (jv Glave’ . y AGRE peG J, ) 
4j=0 ; 1? oe 
Vv me De Gi, )] 
fe phe ne iteo re | (5.14) 


where v(ktm) is the response of the linearized system between sampl- 


inesinetants.. Yor this particular licase it jis eiven by: 
k 
vein) = 0g (etm-j) fy) - vG)] (5.15) 
j=0 
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Equation (5.14) gives the response of the nonlinear discrete- 
data system between sampling instants when the system does not contain 
a ZOH device. 

The series corresponding to (5.14) for the case of systems 
With ZOH is Agee by replacing v(k) in (5.14) by atk) from equation 
(45. 70)8 
De weno Lao Ly Considerations 

Let us consider first the stability of the linearized system. 


The response of this System between sampling instants is given by: 
k 
vikim) =) g (ktm-j) Ly) - y,v(3)1 (5.16) 
j=0 


It will be assumed that the linearized system is stable at sampling 


fistants, chen) stability of (6). 1s" agsuced, if: 


y. |g, (etm) | < © (Say) 
k=0 


which implies that 
lim |g (k+m)] = 0 (5.18) 
O 
ko 
Note that, by stability we imply BIBO stability as well as asymptotic 
stability. 

Now, we proceed to the case of nonlinear systems. Stability 
of nonlinear discrete-data systems follows directly from the stability 
of these systems at sampling instants, and the stability of the linear- 
ized systems between sampling instants. This can be shown as follows: 

The intersampling response of nonlinear discrete-data 


systems is given by equation (5.7), viz.: 
k 
(etm) =)" gi Cktm-j) Ey) - ¥yx(4) — 616609] (5.19) 
j=0 
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Substituting for 04, ) one gets: 
k M 
x(kim) =) go (chm) (yG) - 14x) - yr GI (5.20) 
j=0 i=2 


Let us now test for BIBO stability. According to our assumptions, 


sup |x(k)| < © (321) 
O<k<o 

sup |y(k)| < (5.22) 
O<k<e 


and inequality (5.17) is satisfied. By taking the absolute value of 
Detiessdes or. (3.20) we gers 


|x(k+m)| < { sup ]y(k)| + ly! sup |x(k) | 


O<k<e O<k<o 


M 


+ ly lt sup [x10 ) [a ccm) | 


0 <k<o 
7 < 
- oa (5.23) 


Gleariyeetic ww .l.o. -OLE(),25) 1sebounded, and since, (5.23) 1s true, for 
al Siesard em), . Oma, it tollows that 
np |x (k-+m) | < o (5225) 

O0<k<o 

QO<m<1 
which means that the intersampling response of the nonlinear system 
Te pounded. lt) also means thet the series (5.14) is convergent. 

Similar argument can be used for testing asymptotic stability. 

5.4 Illustrative Example 
In this example we will consider the same system considered 


in the example of the previous chapter. The system is shown in 


Fig. 432. The transfer function.of the forward patheis given by. 


Ts 


om 
G(s) = See RY (5595) 
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The intersampling impulse response of the forward path will be obtain- 


ed in the following way: 
The modified z-transform of G(s) is§ Found’ to be: 


7 (5376) 


z(z- e ) 
Taking the inverse modified z-transform of equation Ce 28)" we get: 


g,(a-lm) = 2-*[G (2,m)] 


i 


-1 
>. Residues of Go (zsm)z" at the poles 


of G (zm)z" > (5.2/7) 
(See reference 38). 
For n=), Go (zm) has two poles at z=0 and eet SO enat: 


og ees —mT ae es 
—T + (l-e )- ar) ee oe 
e e 


Il 
{ 


g, (O+m) 


= (l-e ) (5.26) 


=i 
Fort n>, Go (z,m) has only one pole at z=e , so that: 


-T. -T(n-2) 
e 


g(n-lim) = e (l=en ye (5229) 
Letting k=n-1, we get: 
g (ktm) = pis (loens ec etngecnn skola ee (5.30) 
Now, for T=1 second, we get: 
1- Aa for k=0 
g (ktm) = (So!) 


m 


Ot aevonsesy ace for k>0 


Now, 


y |g, (ktm) | =|1-e "|+ 0.632 sey (GSe3)cme 
k=0 k=1 
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It was shown by equation (4.77) that the linearized system 
is stable at sampling instants. Also equation (5.32) shows that in- 
equality (5.17) is satisfied. Hence, the linearized system is stable 
between sampling instants. Furthermore, the nonlinear system was 
shown to be stable at sampling instants in the region in which in- 
equalttres (4.96) "and (45100)"are satisfied.” "Hence, it»ftoltows that 
the nonlinear system is also stable between sampling instants. 

It has been found in chapter IV that the bound on the in- 
pueey’ (tt), inorder that the nonlinear system or Fie. *452 "is" BiBO 


Stable at sampling instants, is given by inequality (4.106), viz.: 


Vee = | |y*ce) | | = sup ly (k) | << IL, 1k63 Clas = 
O<k<e 





x(0) |) (Seas) 


The corresponding bound on the output of the nonlinear system at 


sampling instants is found to Bae 


Kee eC) |= sectip exe) | Pe 255 | (5.34) 


O<k<a 
A bound on the intersampling response of the driven nonlinear 
system is now readily obtainable from inequality (5.23). Substituting 


(Gro. meron 5704), into (5,23) this bound is) found \to be: 


= sup. |x(ktm)| <[5.435.- 1.163|x(0)|J[|1-e, |te, J (5.35) 
O<k<co 


O<m<1 


x* 
m 
If the system is undriven, then Y*=0 and, provided that (4.107) is 

satisfied, the nonlinear discrete-data system is asymptotically stable 
between sampling instants. The bound on its intersampling response in 


this case is given by: 


x* < 4,06 [la-e “|+e "] (5.36) 
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CHAPTER VI 


CONCLUSION 


In this thesis we have dealt with the response of OGe ee 
discrete-data systems when the initial conditions are not zero. We 
have also pAb the stability problem under these circumstances. 
After reviewing the main results in the field of nonlinear discrete- 
data system stability, non-zero initial conditions are included in a 
general expression for linear discrete-data system response for both 
systems considered. Bounded-input bounded-output stability as well 
as asymptotic stability of linear systems have also been discussed. 
The results are ees to the case of nonlinear discrete-data systems, 
aoe by solving the system response equation by iteration, a discrete 
Volterra type series is obtained in which non-zero initial conditions 
are accounted for. Uniqueness and convergence of the series solution 
obtained in this manner, as well as bounded-input bounded-output 
stability of the nonlinear system, are investigated Sees use of the 
contraction mapping theorem (Banach fixed point theorem) of functional 
analysis. By equating the input of the system to zero, asymptotic 
Stability of ane nonlinear system is established in the region where 
the conditions of the contraction neers theorem are fulfilled. An 
example is given to ine ee the results obtained. 

Finally the response of nonlinear discrete-data systems 
between sampling instants is considered. The response takes the form 


of an infinite series, which is shown to be bounded. Bounded-input 
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bounded-output stability as well as asymptotic stability of nonlinear 
discrete-data systems between sampling instants are also investigated. 
It has been shown that the region of system stability between sampling 
instants is identical to the region of system stability at sampling 
instants, i.e. the bounds on the input and initial conditions are 
identical in both cases. 

It is worth noting here that since the contraction condition 
is a strong requirement, the region obtained by this method are bound 
to be on the conservative side. However, unlike the method used by 
Alper, this method has the advantage that the series obtained is a 
convergent series. 

In conclusion the author wishes to state that the main con- 
EGIDUELONNOTMEN I SP thestoers. Bie kimelusionbor tion=zerominitialecondi= 
tions in the Volterra series representation of the response of a class 
of nonlinear discrete-data systems at sampling instants, and the 


determination of system response between sampling instants. 


ee 








TWorniinon ‘eo ql eae plaodgmya ie Tow o& BRR 
» baa gi eavyol ovle vy init hci galt ene Ta ousiog nade a 
S113 pre 2 WoO » wei | peg Pus My To ny tgs of surly nee, asod 2 


aniiqmme de yolildeiy megeee Go no tea She on Dao kanebd al 
ash ert iewas Leloivl one Seq ata no site all 1.2 snag 
.2o0h2 ftoad nt fear snobe * a 
nudai hoe asters f) 4yite” diet ocd gpesen wow abet 7 : 


Laued sie heron elyo yd bonletdo polhey B19 . See See ee 
1 voy ibdatdam ote sth ban eta aie eure slo etsens ofy @o sf 63 
untae y ie Bot ois gee) ooadaeyls ape wed fowen ebds ,S0TA 
-e3h708 Jaegyeynos y 

iw vedtyys. silo qodentsages af - 

rie ul ates oon J, 90?) aulied od4 2b apie Sie ie eee 

en BOF to Is pe ery (oi tae Aaa a9 Re Bae 


77 ; ') fom love 2obete 360h vebald nen Je 7 


1.0 moyiod wetphor gqegeye Tn wibieiaedes —_— 





60 
REFERENCES 


W. Hahn, "On the application of the method of Liapunov to differ- 
ence equations", (in German), Mathematische Annalen, Vol. 136 


b) 


Pp. eou-441, 1953, 


W. Hahn, "Theory and application of Liapunov's direct method", 
(translated from German), Prentice-Hall, Inc., Englewood Cliffs, 


Midi 1a 


Rin. Kal@eangend.).—. pertram, ‘Controlesystem analysis and desien 
via the second method of Liapunov; II Discrete-time systems", 
Transactions: ASME, Journal of Basic Engineering, Vol. 82, pp. 394- 


400, June 1960. 


Kall0catay #StateYspace analysisiotteontrolvsystems'! } (Prentice=Hall, 


Incweaneticwood Clittie wen ding 19674 


Ya. Z. Tsypkin, “On the global stability of nonlinear automatic 
sampled-data systems", Doklady Akademii Nauk SSSR, Vol. 145, No. 1, 


Dp. ovo, July L962, 


Ya. Z. Tsypkin, ‘Absolute stability of equilibrium positions and 
of responses in nonlinear sampled-data automatic systems’ , 
Aveomatika 1, Telemekhanika, Vol’ 245 No. 12). pp. 1607-1615, 
December 1963. 

Ya. Z. Tsypkin, "Fundamentals of the theory of nonlinear pulse 
control systems", Automatic and Remote Control (Theory), Proceed- 
ings of the Second Congress of the IFAC, 1963, pp. 172-180, 


Butterworth, 1964. 


at 









es 4 
+i DUE TR ‘ ; 
7 oe 
| 7 u . pit i ety lenge ons no” "gine z) 
1 ‘ - 
ened if Dal). eapiien pe f 
rf : Thn-abs 
Te 9 Sa 
§ 
Honrd | Lei 
f j 
{ st * 
4 
* ; 18 
+ 
‘ . 
> atx 
, oe” ; 
' s . 
‘ . 
; 
— 
wl f “7 Yo. a 
os 















é 4 | , » 4 She ¥ _ W : Y T j Ea tj . i detec <iiee 
i . -EOCL ve San 1390 
ee oT = : 
-_ pélog YBeis fries ek (pt) 203520 elastin nhunwa” ieee, 
id a 





ame as , a rp, tam 9. nh" e win : re an ce) 
ais ‘i on ie jane) dares: a . . 


LO® 


1 


IDR 


Ss 


EA. 


15. 


L6. 


61 


Ya. Z. Tsypkin, "Frequency criteria for the absolute stability of 
nonlinear sampled-data systems", Avtomatika i Telemekhanika, 
Voitee) .eNochs.Pppre234-209-iMareh 1964. 

ExieJuryoend B-W. Lee, “On the stability of a certain class of 
nonlinear sampled-data systems'', IEEE Transactions on Automatic 
Comer olymVolseAGz 9 eNo wale po, 51-61.) danuary,1964% 

Ee.is Jury and’B.W. Lee,, "On: the absolttevstability of nonlinear 
sampled-data systems", IEEE Transactions on Automatic Control, 
Vor, “AG-94eNo. “Sa nnse551—-554en0ctobers1 964% 

Beds@Jury andsBeW. Leese Thefabsolute stabilitzy ofasystemsewith 
many nonlinearities", Avtomatika i Telemekhanika, Vol. 26, No. 6, 
pp. 945-965, June 1965. 

E.I. Jury and BAW. Lee, “A stability theory for multinonlinear 
control systems", Proceedings of the Third Congress of the IFAC, 
1966 .fopte 2eAnd=28A0 U1, 

G.W. Johnson, "On Liapunov stability vs. bounded input~bounded 
OUP LEER Transactions: on Automatic Control, Vol. AG=9, No. 25 
pO mivo-/9% Apraigloo4. 

L.A. Zadeh and C.A. Desoer, "Linear system theory", McGraw-Hill 
BookiConss NewivoxckSONey., 1963. 

I.W. Sandberg, "Some stability results related to those of 

ViM. Popov'';. Bell Se eeen Technical Journal, Vola 642eppesZisee 
2148, November 1965. 

A.R. Bergen, R.P. Iwens and A.J. Rault, "On input-output stability 
of nonlinear feedback systems", [EEE Transactions on Automatic 


Control; Neon, Nowy poee/42-744, October 1966. 


iA 





' 3o eenlo pleogss s 30 wiEdade ante-q0" (ond it eal an 
StiomeguA oo voels comet? Gaal y Bitte ee nJalhet quiba yeonl non 
pol yoauent , fede .oq .f -o8- O60 ie fee 
Yhonifnen. 7H yrl ble te sbhineds no) 00° ,&tl «Mee Be gol Ie Or 7 
4092069 zs) suwotul ne eyo iaebeiey 7 , emeseve Pere ee Aa 
OL sgdodoanh , eo 182 .9g 8 20 Sie . ; , 
. on i Jiev aa wioT F ativtougeh ,"eals heenhiaod Saee 


. 
fjiv poadewe Jy qilltdotke sauloedaé oft” .set . Well Gp We on ee eee 
B00) omit. ,2a-ANe gq : 
ive tilidsss A" aod Wel hrs Gl eT aE 
) Lib! ovo 7) sgnifessev? ,"amedese Tesahes 
LL ABS LARS woe , dO 
SPs curs hebdne Vivi bids de eas ae id pO" , noetol We, EE 
| aiohassaueyT SW0t , “apqgeo 


PORT LPgA pOUE-BUL cog 
repay » Vibes th tbo! ecen@ -AVO chan Mebat Mae aL 


-COPL orl ih gato’ Wa" 5.00 Bead 
" . <4 . 
iu SP is 7) ree | r, ail ford Wd Dee ie pith readin «Wal. tL 7 
ie 
SELES agg Ad slov. .luntwiatk footaiset matey Lhe ,"Vuged a0 ae 
: ; es ‘ 





Lie 


ioe 


19. 


DOR 


2 


Jee 


Zoe 


24. 


2S. 


62 


R.P. Iwens and A.R. Bergen, "On bounded-input bounded-output 
stability of a certain class of nonlinear sampled-data systems", 
Joumnalwot the erankiin institute, Vol. 282, No.4, pp. 193-205, 
October 1966. 

R.P. Iwens and A.R. Bergen, "Frequency criteria for bounded- 
input bounded-output stability of nonlinear sampled-data systems", 
IuEE@lransectionsiion Automatic Controls) Vol¥ AC2I28.No. 1, 

Pp. £0 —ewer ebruary .196/). 

JgbarBarrett, ‘The ase of Wolittersatseries ito find region! of 
stability of a nonlinear differential equation", International 
Joma womucontrol, Vol. 1. No. 3, pp. 209-216, March 1965. 

G.S. Christensen, “On the convergence of Volterra series", IEEE 
Teansact tons zon’ AutomatvenGontroly Wolds” ACHioysNoeG Epp. 736-737, 
December 1968. 

G35. 4 Ciigistensenmand Gewetrott,, (Onethe Wnclus ion ior fad tial 
condieons! in Volterraisertes!.. International ournal cf Control, 
In press. | 

P. Alper, "A consideration of the discrete Volterra series", 

IEEE Transactions on Automat ite) Contnolevo AARC, Nomis. 

Di Boe 2e 32/5 Jul 96 5e 

R.S. Rao, "Stability studies on nonlinear discrete-data systems 
through functiondlsanalysis!,) Ph. DiMThesis, University of Alberta, 
Edmonton, Alberta, Canada, Fall 1969. 

ReSeukao and G.SsuChristensen, 7.00) the convergence of wdiscrete 
Volterra series", IEEE Transactions on Automatic Control, Vol. 
AG=1 5), “Nose 14) February 97/0). 


T.L. Saaty and J. Bram, "Nonlinear mathematics", McGraw-Hill 


Book Co., New York, N.Y., 1964. 


Sa 







9124 Beicre bot tppal Seamti-babaond 00” mgr! yA aes 
‘aieieye wisteintenan shoal ion Je Wiels aber te waehitdnses 
Hef eP|! da.e _cb2 2187 ,agesisenl wiisga wifd! Ye: Inpisct 

\ | O901 tedatad: 
bhbived vot slits Yenaypes 9" , ose \ aye ee See 
. fqoag thet load So yaelicete tne Gal i 
1... CuY , Tavined olsewadah ag emigedanayT SSN: 
‘ol eectae g6t-o8 s9@ 
VY io »al ef  ~agepse® 00° lee 
PIsaes fh aiottbionn » Jo wehT dese 
lus “Jf | [Ov hee eew®> Se. Daneel 
" wensigieds .23 JOS 
; é Ve Mee tik eno kive anne. 
Meet Tredmeoed 
. heat -cparstelrd. 12.9 8 


pyisaieY ak emls)baod 


-ea5%q al | 
; ‘yy Si uo 4] lat ot eerie a” , 2S0lA 7! $k 
POP Wh aeal Leh (ti) ewok Jon eeB eT taal 


20D vfot  <Le-o8e “7a ; 

RMAs S3hleal oe iagkl tear go ealinces Yahi” pee WALES ee 
(OTe 29 eka peluedl 20M» ebae tam Lennd toad? gaokda.. _ 
; OY. L1H, | | 


: oy a yi O ae id 7 a ? a A be f 





; ‘ 7 
7 on B aes — ; : 
. 7 
. 


63 


26. L. Liusternik and V. Sobolev, "Elements of functional analysis", 
(translated from Russian), Frederick Une aeePubl shine “Cox,, 
NewiYouk,, Nw, 96d. 

2H E.T. Copson, "Metric spaces", Cambridge University Press, London, 
England, 1968. 

2c. J.€. Hsu and A.U. Meyer, "Modern control principles and applica- 
tions", McGraw-Hill Book Co., New York, N.Y., 1968. 

29. G.P) Szego and J.B. Pearson, Jr., “On the absolute stability of 
sampled-data systems: the indirect control case", IEEE Trans- 
aeeLOusmon Automatic Control, Vol. AC-9, No. 2, pp. 160-163, 
April 1964. 

Spe oeeepel'man, ‘A Criterion for the absolute stability of non- 
linear sampled-data automatic systems", Avtomatika i 
Tewemekhanikva,. VOlwn 2/7, No. 10, pp. /4-/7,,. October 1966. 

Ble Ms DeRusso,—R.J. Roy and C.-M. Close, "State variables for 
eueineers. = Jonn Wiley & Sons, Inc., New York, N.Y., 1967. 

a2 Vaeuz- [sypkin, Sampling systems theory and its applications. , 
(translated trom Russian), Ihe MacMilian Cow, New, York, 1964. 

33.  A.E. Taylor, “Introduction to functional analysis", John Wiley 
SEOOts lnc.) News YOrk waNeis ,. 1956. 

BAe W.A. Porter, "Modern foundations of systems engineering", The 
Machiilan Co., New York; N.Y., 1960. 

555 E.I. Jury, "Theory and application of the z-transform method", 
John Wiley & Sons, Inc., Near NMordles Wipes ey TEU 

B65 E.l. Jury and F.J. Mullin, "A mote, on the operational solution 
of linear difference equations", Journal of the Franklin Institute, 


Vol. 266, No.93, pp. 69-205, September 1958. 


tA 


: i 
, Blevioun Jnavitom’d ia gineoolW" yvhhadad LV bie nena i 
‘i 
, ©” galvalldull 46g) anivehor? | Gnteep eee beouluants) 
fant cet Prt elt 


ei hot, | ¢# ari § iwoasV ir ost invita 4 sa a oa atyool” histo. ts 
POO). ~hoal gent 






Hos T yal ) 1 a4 Fala | la Te mah” eye 0 ih keeate arall Wek 
Yo’ ,aoY wall ..99 due! (19 iegpabaily “enka 


. . “4 Th  coemeel 0) Lie Geee® 187 ee 


ipo» 2oett bah af) alee epee eheleua 

{Ty ¢ lop ,e=)) rev, mo? Si jp wo enpigos 7 
-O0@L Livga | 

71 insti A". nna feet 2a 88 


Kats , ve oi levoe @ehiieleug@ tegasf°~ 


loVY ,etinedtonsl or 
LR, oeawied HeG IE 


Jiv adel ,"wementeae 


nity?” oP aT: f ae¥ SE 
(meee gia he?elews12) 
_ rik { ot oe: Geis busin”. ander .oe €E 7 
' ) Ofel ,.0.4 . 266) tek sank panel a : 


ad T ‘ ; fia “4250 A? { ee el 4 y 9 } inf hy P| iain | i val" ' 149304 7 i he 


O0RL:,.Yot Set aah eles wih tin hadi - 


BONIS NiolSiis4~. 943 Yo maha leing nn ee 
yt set a Wuti¥: . L 





Bes 


De 


Ue 


40. 


64 


MiB, Brilliant, “Theomyrontihe analysis of nonlinear systems", 
MEL Research Laboratory of Electronics, Report 345, 1959. 

B.C. Kuo, “Analysis “and synthesis ‘of sampled-data control 
Systems , Prentice-Hall, Inec.; Englewood Cliffs, N.J. 1963. 
PAGaeDOLe. time domain analysis and design of control systems", 
Addison-Wesley Publishing Co., Inc., Reading, Mass., 1965. 

J.T. Tou, "Digital and sampled-data control systems", McGraw- 


Pie eb OO Ke CO, eNews VOLK aN ai sts LI OO 


aa 





e"sansey? wcclined Je @leetene of9 46 vient" 4 
UC0F ,2eo syveqal ,exdaevitdt? to oti ieee Te 
bovine atnh-bsiqmae Yo ehialieg tne ehaytnoa* ot 0.8 

| 000 VL eID boototget.,. 27 , Uel-aghgass§ eeesave 
,"amereye loro Ww ngtesh beg efovinen winasb oabt™ ,270% .3s2 

‘COLD , exalt quehlia sh ae »-09 pol dalide® valenenes 1 BGA 

watGut] ,"omadeve losnnwe efeb-hateane bus fastera" wot ‘TB 


Pre Vo deet eek et doad 1 





65 


APPENDIX 


Z-TRANSFORM THEORY 


The theory of the Z-transform is well known. Its develop- 
ment from the Laplace transform and its main theorems are discussed 


in the textbooks on the subject of sampled-data Seateng a How- 


ever, the defining equation of the Z-transform will be derived here, 
then its main theorems and the different methods for the determina- 
tion of its inverse are stated. 
A.1l The Z-Transform 

A uniformly sampled signal x*(t) of uniform sampling period 


Py ts’ a *sienal “such ‘that 


as) . § (een) 
n=0 


co 


eS Cnt) *OCrenD) | Ge) 
n=0 


I] 


x*(t) 


where x(t) is the continuous signal before sampling. This means that 
a sampled signal is a train of impulses separated from each other by 
T seconds and each is of height equal to the value of the sampled 
signal at the sampling instants. It is assumed that the continuous 
Sion eect) aise zero Hor *t-0% 


Tf mow equation (A/1) is Laplace transformed, one gets: 


ioe) 


x*(s) = ye (al) Borne (Az?) 
n=0 


Ts : : 
Letting 6 “=z, é6quation “(A.2) may be written as: 
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oo 


x*(s) l = x2) 2» <(n1) 3 (A. 3) 


s=—l]nz 
Al n=0 


Equation (A.3) defines the Z-transform of the signal x(t). 


A.2 Theorems) of the Z-Transform 


In this section only the theorems used in the development 


presented will be stated. 
A.s2.1. Linearity 


ies 8, (t) and 8, (t) are Z-transformable functions and if 


a and b are two constants, then 


Z2[ag, (t) + bg, (t) ] =a Z{g,(t)] 2b Z{g,(t)] (A. 4) 


where 2Z[ ] represents the process of Z-transformation. 
A.2.2 Real Hansdatien 

If g(t) is Z-transformable function and if n is a positive 
integer, T is the sampling period, then 


Z{g(t -nT)] = 2 Z[g(t)] | (A.5) 


and n-1 


Fie se) me wee wie oa] y e(eiyen | (A.6) 
k=0 


A.3 The Inverse Z-Transform 

The determination of the sampled signal from its Z-transform 
is as important as the determination of the Z-transform of a sampled 
signal. This process of passing from the Z-domain into the time 
domain is called "the inverse Z-transformation". | 

To evaluate the inverse ieee one of the following 


three methods can be used. 
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Ae eee ine inversion Rone lat 


The sampled signal g*(t). defined by 


ee) 


gk(t) = g(nT) 6(t-nT) (A.7) 
n=0 , 


T being the sampling period, can be evaluated from its Z-transform 


G(z) by the use of the inversion integral 


aCe ae p CC) meena (A.8) 


i 
where [ is a circle of radius age centered at the origin of the 
z-plane, and c is a constant such that all the poles of G(z) are 
enclosed by the .circle. 
By the use of Cauchy's residues theorem equation (A.8) can 
be rewritten as: 


g(nT) = Sum of the residues of G(s) hes 


ate the=poles: of, G(z) a (A.9) 


A.3.2 Power Series Expansion 


Here the Z—transform G(z) is written in the form: 











(A.10) 


where for physical systems m <n.’ 
By performing long division G(z) may be written in the form: 


al 2 sah 
= 4 ooo 4 4s ap eee ACLL 
Ce) ana 2 2a ot + a2 ( ) 


Now by taking the z-transform of equation (A.7) term by term one gets: 


Cea) 2 MG) GD 2 OU een een eee Gn) 


* see in particular pp. 66-68. See also Ref. 40, pp. 176-179. 
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By comparing equations (A.11) and (A.12) we can conclude that when 
G(z) is expanded into a power series of sed the coefficient of 
z in the expansion corresponds to the value of the signal g(t) 


at t=nT. 


Awe oeeeartial Fraction Expansion 


Gz) 
Z 





Here is expanded into partial Fractions of the form 





where the a,'s are constants and the c,'s are the poles of had 
ue aya 
Then G(z) will be formed of terms of the form 





By the use of 
i 
z-transform tables the time function corresponding to each term is 
Obtained to get a time function g(t) corresponding to G(z). This 
time function, however, is of significance only at sampling instants. 
It should be noted that any of the above three methods 


evaluates the sampled signal g*(t). They do not give any information 


about the behaviour of the signal g(t) between sampling instants. 
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